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Abstract 

We present a calculation of the matrix elements of the most general set of AS = 2 
dimension-six four-fermion operators. The values of the matrix elements are given in terms 
of the corresponding 5-parameters. Our results can be used in many phenomenological 
applications, since the operators considered here give important contributions to K^~K^ 
mixing in several extensions of the Standard Model (supersymmetry, left-right symmetric 
models, multi-Higgs models etc.). The determination of the matrix elements improves the 
accuracy of the phenomenological analyses intended to put bounds on basic parameters of the 
different models, as for example the pattern of the sfermion mass matrices. The calculation 
has been performed on the lattice, using the tree-level improved Clover action at two different 
values of the strong coupling constant {(5 = 6/go(a) = 6.0 and 6.2, corresponding to o^^ = 2.1 
and 2.7 GeV respectively), in the quenched approximation. The renormalization constants 
and mixing coefficients of the lattice operators have been obtained non-perturbatively. 



1 Introduction 



Important information on the physics beyond the Standard Model (SM), such as supersym- 
metry, left-right symmetric models, multi-Higgs models etc., can be obtained by studying 
FCNC processes. Among these, AF = 2 transitions play a very important role. They have 
been used in ref. ||^, ^], for example, to put constraints on the sfermion mass matrix. In this 
paper we present the results of a lattice calculation of the matrix elements of the most gen- 
eral set of AS* = 2 dimension-six four-fermion operators, renormalized non-perturbatively 
in the RI (MOM) scheme [^-[0]. Our results can be combined with the recent two-loop 
calculation of the anomalous dimension matrix in the same renormalization scheme [^] to 
obtain K^-K^ mixing amplitudes which are consistently computed at the next-to-leading 
order. A phenomenological application of the results for the matrix elements given below, 
combined with a complete next-to-leading order (NLO) evolution of the Wilson coefficients, 
will be presented elsewhere 1^]. 

mixing induces the neutral kaon mass difference AMx and is related to the 
indirect CP violation parameter ek- In the Standard Model, this transition occurs via 
the dimension-six four-fermion operator O^'^^^, with a "left-left" chiral structure. The 
-B-parameter of the matrix element {K^\0^^^'^\K^) ^ commonly known as Bk, has been 
extensively studied on the lattice due to its phenomenological relevance p^] , and used in many 



phenomenological studies |11]. For the other operators, instead, all the phenomenological 
analyses beyond the SM have used S-parameters equal to one, which in some cases, as 
will be shown below, is a very crude approximation. The present work is the first fully 
non-perturbative study of the matrix elements of the complete set of AS* = 2 four-fermion 
operators. With respect to other calculations, the systematic errors in our results are reduced 
in two ways: i) by using the tree-level improved Clover action [^] and operators [^] we 
obtain matrix elements for which discretization errors are of ©(a^a); ii) by renormalizing 
non-perturbatively the lattice operators, we have eliminated the systematic error due to the 
bad behaviour of lattice perturbation theory. With the non-perturbative renormalization, 
the residual error is that due to the truncation of the continuum perturbative series in the 
evaluation of the Wilson coefficients of the effective Hamiltonian. This error is of 0{a1), 
where a<j is the continuum renormalization constant and could in principle be reduced by 
the calculation of the N^LO corrections in the continuum. In the RI scheme, for example, 
this step has recently been done for the calculation of the renormalized quark mass jl4{ |. 

For the reader who is not interested in technical details, we now give the results for the 
-B-parameters of the relevant operators. The choice of the basis is arbitrary, and different 
bases can be found in the literature, see for example |||] and @. We have used the SUSY 
basis which is also the one for which the numerical values of the Wilson coefficients, computed 
at the NLO, will be given Q, namely ^ 

Oi = s"7a.(1-75K AA.(l-75)d^ 

03 = s-"(l-75)'i''^^(l-75K, (1) 

04 = r(l-75)d" s^(l+75)d^, 



^ We use here the Euclidean notation. 
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where a and /? are colour indices. The i?-parameters for these operators are defined as 
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(2) 



where the notation Oi{fj,) (or simply Oj) denotes the operators renormalized at the scale 
/i. A few words of explanation are necessary at this point. In eq. (^) operators and quark 
masses are renormalized at the scale in the same scheme (e.g. RI, MS, etc.). The numerical 
results for the S-parameters, Bi{iJ,) computed in this paper refer to the RI scheme. Moreover, 
without loss of generality, we have omitted terms, present in the usual definition of the Bs, 
which are of higher order in the chiral expansion. Since the definition of the i?-parameters 
is conventional, we prefer to use those in eq. (Q) for which, as explained in sec. ^, the scaling 
properties are the simplest ones. 

Our best estimates of the S-parameters, for a renormalization scale of ^ = 2 GeV are 



53(/i) 

B^iix) 



0.69 ±0.21, 
0.66 ±0.04, 
1.05 ±0.12, 
1.03 ±0.06, 
0.73 ±0.10 . 



(3) 



The remainder of the paper is organized as follows: in sec. |2|, we address the problem 
of operator mixing and renormalization and give a brief account of the non-perturbative 
method (NPM) used in the computation of the operator renormalization constants; in sec. ^ 
we discuss the definition of the i?-parameters and describe their extraction from the lattice 
correlation functions; in sec. ^, we present our results for the full operator basis, sX [5 = 6.0 
and 6.2 and for different renormalization scales; a discussion of the errors assigned to the 
final results in eq. ^ can also be found in this section; finally, in sec. [s], we present our 
conclusions. 



2 Non-Perturbative renormalization 

In this section, we briefiy recall the reasons for which the non-perturbative renormalization of 
the lattice operators is important and describe the procedure which has been used to obtain, 
for the cases of interest, finite matrix elements from the bare lattice operators. 
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The Wilson lattice regularization breaks chiral symmetry. This implies that each operator 
in the AS" = 2 Hamiltonian mixes with operators belonging to different chiral representa- 



tions 1 15, 16]. Because of the mixing induced by the lattice, the correct chiral behaviour 
of the operators is achieved with Wilson fermions only in the continuum limit. This rep- 
resented a long-standing problem in the evaluation of Bk |17, |18| 0, only recently solved 



with the introduction of Non-Perturbative Renormalization methods. In these approaches 
the renormalization constants (mixing matrix) are computed non-perturbatively on the lat- 
tice either by projecting on external quark and gluon states (NPM) as proposed in ref. |Q 



or, in the spirit of ref. ||16||, by using chiral Ward Identities Recent studies of the 

-B-parameters, with both non-perturbative renormalization methods, and [^l|, show 

that discretization effects are less important than those due to the perturbative evaluation 
of the mixing coefficients. 

Given the success of the non-perturbative methods in the computation of Bk, the NPM 
has been applied to the evaluation of the two AI = 3/2 i?-parameters of the electro-penguin 
operators, By and Bg (these S-parameters coincide with those of the operators O4 nd 
O5 respectively). Also in this case, as shown in ref. Q, it has been found that the non- 
perturbative renormalization of the lattice operators gives i?-parameters that significantly 
differ from those renormalized perturbatively |6|, |2^ . 

An extensive study of the renormalization properties of the four-fermion operators can be 
found in . There we detail the issues of relevance to the non-perturbative renormalization 
of all the A^* = 2 operators. We have used these results in the present study. 

The NPM for the evaluation of the renormalization constants of lattice operators con- 
sists in imposing suitable renormalization conditions on lattice amputated quark correlation 
functions In our case, we compute four-fermion Green functions in the Landau gauge. 
All external quark lines are at equal momentum p. After amputating and projecting these 
correlation functions (see refs. @ and for details), the renormalization conditions are im- 
posed in the deep Euclidean region at the scale p"^ = ^x^. This renormalization scheme has 
been recently called the Regularization Independent (RI) scheme [p3|| (MOM in the early 
literature) in order to emphasize that the renormalization conditions are independent of the 
regularization scheme, although they depend on the external states used in the renormal- 
ization procedure (and on the gauge). Thus, at fixed cutoff (i.e. fixed /3), we compute 
non-perturbatively the renormalization constants and the renormalized operator O^^ (fi) in 
the RI scheme. In order to obtain the physical amplitudes, which are renormalization group 
invariant and scheme independent, the renormalized matrix elements must subsequently be 
combined with the corresponding Wilson coefficients of the effective Hamiltonian. For the 
operators of interest, the latter are known at the NLO in continuum perturbation theory |^]. 

In 0], we have determined non-perturbatively the operator mixing for the complete basis 
of four-fermion operators, with the aid of the discrete symmetries (parity, charge conjugation 



and switching of flavours), in the spirit of ref. [17|. The renormalization of the parity-even 



operators, relevant to this work, is more conveniently expressed in terms of the following 
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In the staggered fermion approach, where chiral symmetry is partially preserved, the AS ~ 2 matrix element 
displays the correct chiral behaviour. Thus, the i?A'-parameter obtained with staggered fermions | [l9t has been 
deemed more reliable. 
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basis of five operators: 



= V + A-x A, 
Q2 = V xV - Ax A, 

Q3 = SxS-PxP, (4) 
Q4 = S xS + PxP, 
Q5 = TxT. 

The operators Qi, . . . , form a complete basis on the lattice. In these expressions, TxT 
(with r = V, A, S, P,T a generic Dirac matrix) stands for ^(V'ir^2'03rV'4 + V'irV'4V'3rV'2); 
where tpi, i = 1, . . . , 4 are fermion fields with flavours chosen so as to reproduce the desired 
operators (see ref. for details). 

The parity-even parts of the five SUSY operators defined in eq. (g), which are the relevant 
ones for K^-K^ mixing, are related to the operators of eq. (Q) in the following way: 

01 = Qi, 

02 = Qa, 

03 = -^(Q4-Q5), (5) 

04 = Q3, 

05 = -^Q2. 

On the lattice, Qi mixes under renormalization with the other four operators as follows 



Qi — Zu 



5 

1=2 



(6) 



where Zn is a multiplicative logarithmically divergent renormalization constant; it depends 
on the coupling and afi. The mixing coefficients Zu (with i = 2, . . . , 5) are finite; they only 
depend on the lattice coupling gQ{a). 

The other renormalized operators are defined as follows: 

Q2 = Z22Q2 + Z23Q3, 

Qs = Z32QI + Z33QI (7) 

Qa = ZaaQa + Za^QI-, 

Q5 = Z54QI + Zr,^Ql, 

where the ZijS are logarithmically divergent renormalization constants which depend on the 
coupling and a/j,. The above mixing matrices are not peculiar to the lattice regularization, 
but also occur in the continuum. The breaking of chiral symmetry by the Wilson action 
requires the additional subtractions: 

Qi = Qi + ^ Z^jQj, i = 2,3, 

Qt = Qi+ Yl ^tjQj, i = 4,5, 

i=l,2,3 
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where the ZfjS are finite coefficients which only depend on gQ{a). 

The results for all the renormalization constants Zij and Zfj (computed with the NPM, 
for several renormalization scales at P = 6.0 and 6.2) can be found in [^. 

3 B-parameters 

In this section we discuss the definition of the i?-parameters and their dependence on the 
renormalization scale. We also sketch the extraction of these quantities from lattice correla- 
tion functions. 

The i?-parameters are usually defined as 

, («°|6.(M)I«°> _ 
{K»\Oi\R-')vsA 

where the operator matrix elements in the Vacuum Saturation Approximation (VSA) are 
given by 

{K'\^^\K')vsA = 2(l + ^j\{K'\AM\ 
{K'\d2\K')vsA = -2 (l - \{K'\Pm\ 

{K'\d,\K')vsA = (l- J-)Ki^°imP, (9) 
{K'\6,\K^)vsA = 2|(i?0|P|0)p + ^|(i?0|i^|0)|2 , 
{K^\d,\K^)vsA = ^\{K^\P\0)\'' + \{K^\AM' ■ 

and P are the renormalized axial current and pseudoscalar densities, = Z^A^ and 
P = ZpP, with Za the (finite) renormalization constant of the lattice axial current, = 
s7/i75(i, and Zp the renormalization constant of the lattice pseudoscalar density, P = s'y^d. 
For simplicity, P is renormalized at the same scale and in the same renormalization scheme 
as the four-fermion operators (the RI scheme in our case). Using the relations 

|(/?"|i,.|0)|= = Mj/i, 

where the second equality is a consequence of the Ward identity for the axial current (with 
ms{lj) and m(i{^) renormalized in the same scheme and at the same scale as P), we find, 
with A^c = 3 

{K'^\^,{^i)\K^)ysA = iMlfl, 
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2 f2 



M'kI 



Mlfl 



The VSA values of the matrix elements of O4 and O5 in eq. ( pA| ) differ from the factors 
appearing in the definition of the B-parameters in eq. @) by the terms proportional to 1/6 
and 3/2 respectively. These terms, which originate from the squared matrix elements of the 
axial current in eq. (|9|), are of higher order in the chiral expansion and have been dropped 
in our definition of the i?-parameters. This implies that, out of the chiral limit, the values 



of i?4 and with our definition differ from those obtained by using eq. (11). To illustrate 
this point, let us imagine that as rag 00, for some value of the renormalization scale fi, 
the values of the matrix elements of O4 and O5 were exactly those of the VSA. Under these 
hypotheses, using eq. @, we would get 84(12) = 7/6 and B^{jl) = 5/2 instead of one. 

We now explain why we prefer the definition of the i3-parameters given in eq. (|2|) rather 
than the standard definition obtained from eq. (|TT|). Neglecting discretization errors, the 



i?-parameters of the operators defined in eq. (Q) obey the renormalization group 

equation 



(70. - 2jp)B,if,) 



(12) 



where fj,d/dfi = fid/d/j, + (3{as)d/das, and 70, (tOjOj) and 7p are the anomalous dimen- 
sion (matrix) of the operator Oi{^) and of the scalar density respectively |^. The physical 
amplitude is given by 



{K^\neff\K^) 



Ci{Mw/fi){K''\6i{fi)\K^) 
Ci{Mw/f^) X Bii^i) X 
as{Mw)\-^o,/2l3o 



1 



msip,) + md(/i))2 

(70,-27p)/2/3o 



(13) 



where, in the last expression, we have only shown the leading behaviour of the different 
factors which depend on ^, namely the Wilson coefficient, the i?-parameter and the quark 
masses. Eq. (p^ ) shows explicitly the cancellation of the //-dependent terms in the ampli- 
tude: the quark masses scale with an anomalous dimension which is opposite in sign to 
that of the pseudoscalar density (since m(//)P(/i) is renormalization group invariant) so that 
Bi{pi) /m?{^) scales as the corresponding operator Oi{fi); the //-dependence of the latter is 
then cancelled by that of the corresponding Wilson coefficient. This remains true at all order 
in as- Out of the chiral limit, with the standard definition of the i?-parameters obtained by 
using the VSA matrix elements of eq. (lTl|), the scaling properties of 84(^1) and B^{ijl) would 
have been much more complicated. The reason is that, in these cases, the two contributions 



For simplicity we ignore the mixing of the operators O2-O3 and O4-O5. 
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on the right hand side have a piece which scales as the squared pseudoscalar density and 
another (proportional to the physical quantity |(^'^|A^|0)|^) which is renormalization group 
invariant. The ^u-independence of the final result would then have been recovered then in 
a very intricate way. Since the definition of the i?-parameters is conventional, we prefer to 
use that of eq. (|2|), for which the scaling properties of all the i?-parameters are the simplest 
ones. Moreover, with this choice, they are the same as those derived in the chiral limit. 

In order to extract the i?-parameters, we need to compute the following two- and three- 
point correlation functions: 

Gp{t,,p} = J2{Pix)P\0))e-P^-^ , GAit,,p) = Y.{Ao{x)PHo))e~P-^ (14) 

X X 

G^(tx,ty;p,q) = ^(Pt(y)0(o)pt(^))e-pyv-; 

where x = {x,tx),y = {y,ty) and O stands for any renormalized four-fermion operator of 
interest. All correlation functions have been evaluated with degenerate quark masses. By 
forming suitable ratios of the above correlations, and looking at their asymptotic behaviour 
at large time separations, we can isolate the desired matrix elements 



i?2 



G^^ {K^{q)\d,\K\p)) 



Z\GpGp Z2|(0|P|K0)|2 ' 

1 1 kK\q)\6^\K\^) 

2(1 - GpGp ^ 2(1-2^) zimpm? ' 



{l-^)ZlGpGp (1-J-) ZimP\K^)? ' ^ ^ 



i?4 

R^ 



1 gp, i{K^m6,\K^{p)) 



2ZlGpGp 2 Z|,|(0|P|ifO)|2 ' 



2ZlGpGp 2 Z|,|(0IPIK0)|2 



We stress that the P-parameters extracted from Ri , P4 and R^ are identical to the B- 
parameters for the operators O^'^^^^q^/^ Oj^"^ respectively. In ref. 1^, the results 
referred to the operators 0^^=2,08^^ and Oj^'^ aX (3 = 6.0 only. In this paper, we present 
the results for all the P-parameters and for (3 = 6.0 and 6.2. 



4 Numerical results 

Our simulations have been performed at /3 = 6.0 and 6.2 with the tree-level Clover action, 
for several values of the quark masses (corresponding to the values of the hopping parameter 
k given in table |l|), in the quenched approximation. The physical volume is approximatively 
the same on the two lattices. A summary of the main parameters is given in the same 
table. "Time Intervals" denote the range in time (in lattice units) on which the two-point 
correlation functions have been fitted to extract the meson masses and the matrix elements 
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of An and P. The ratios Ri, related to the matrix elements of the four-fermion operators, 
have been extracted on the same time intervals. Statistical errors have been estimated with 
the jacknife method, by decimating 10 configurations at a time. 

As discussed in sec. |2|, the renormalization constants have been obtained from the quark 
correlation functions, in the Landau gauge. The results for the Zs have been obtained 
on a 16^ X 32 (16^ x 32) lattice at /? = 6.0 (6.2), using a statistical sample of 100 (180) 
configurations. In constructing the renormalized operators we have used the central values of 
the renormalization constants neglecting their statistical errors. For this reason the errors on 
the 5-parameters only include those of operator matrix elements. In the ratios ([l5|), we also 
need the axial-current renormalization constant Za and the /U-dependent renormalization 
constant Zp of the pseudoscalar density. Although Za should not depend on a^u, slight 
variations of its NPM estimate, arising from systematic effects, partially cancel analogous 
variations of Ri, giving more stable results in the extraction of the matrix elements. The 
NPM estimates for Zp and Za used in the present work are those of ref. [p^j. 





Run A 


Run B 




6.0 


6.2 


No. Confs 


460 


200 


Volume 


18^ X 64 


24^ X 64 




0.1425 


0.14144 


k 


0.1432 


0.14184 




0.1440 


0.14224 






0.14264 


Time Intervals 


10-22 


14-26 




42-54 


38-50 


a-^K*) (GeV) 


2.12(4) 


2.7(1) 



Table 1: Summary of the parameters of the runs at (3 = 6.0 (run A) and (3 = 6.2 (run B). The 
calibration of the lattice spacing has been done using the lattice-plane method of ref. W^] . 



In order to extract the S-parameters from the ratios of eqs. (|15|), we follow the procedure 
of ref. by fitting the Ri linearly in X and Y (i.e. linearly in m\ and (p • q) ) with the 
function 

Ri = ai + piX + -fiY , 



where 



with 



SGpGp 3Z^|(0|P|K0)|2 ' 



Y 



(p ■ q) 



X 



m 



K 



p ■ q = E{p)E{q) — p ■ q, sinh 



2fm\ 



V 2 y 

Assuming ai and /3i to be zero, Bi is given by: 

-Bi = 71 . 



sinh 



i=l,3 



2 Pi 
sm (- 



(16) 

(17) 
(18) 
(19) 
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Bi 


B2 


i?3 


^4 


i?5 


0.31 


0.72 ±0.14 


0.70 ± 0.08 


1.62 ±0.28 


1.21 ±0.08 


1.70 ± 0.24 


0.62 


0.70 ±0.15 


0.64 ±0.04 


1.21 ±0.12 


1.08 ±0.05 


0.81 ±0.10 


0.96 


0.70 ±0.15 


0.61 ±0.03 


1.10 ±0.08 


1.04 ±0.04 


0.68 ± 0.07 


1.27 


0.70 ±0.15 


0.60 ±0.03 


1.03 ±0.07 


1.01 ±0.04 


0.69 ± 0.06 


1.38 


0.70 ±0.15 


0.59 ±0.03 


1.03 ±0.07 


1.01 ±0.04 


0.70 ± 0.05 


1.85 


0.71 ±0.15 


0.59 ±0.03 


0.98 ± 0.06 


0.99 ± 0.04 


0.68 ± 0.05 


2.46 


0.72 ±0.15 


0.57 ±0.03 


0.96 ± 0.05 


1.00 ±0.04 


0.69 ± 0.04 


4.00 


0.74 ±0.15 


0.55 ±0.02 


0.87 ±0.04 


0.99 ± 0.03 


0.73 ± 0.04 



Table 2: Values of the B-parameters in the chiral limit, at different scales fx a for (3 = 6.0 



Since we are working in the linear approximation in X and Y, there is no difference between 
the value obtained in the chiral limit and at the physical kaon mass. In the chiral limit, the 
-B-parameters for the other four operators are given by 



(i = 2,3,4,5) 



(20) 



3/2 3/2 

We stress again that B4^ = and = Bj .At the physical kaon mass we have instead 



Bi = ai + {f3^ + -fi)Xs 



(i = 2,3,4,5), 



f21) 



where Xg is obtained by extrapolating linearly X as a function of the squared pseudoscalar 
meson mass to the physical value rn^^ . 

In tables § and ^ we give the values of the i?-parameters in the chiral limit, extracted 
using eqs. (19) and (|^), at /? = 6.0 and 6.2 respectively; in tables ^ and ^ the i?-parameters 
are evaluated at the physical kaon mass using eq. (^l]) for i = 2, 3, 4, 5. At /? = 6.0, the 
results for Bi, B4 and B^ extrapolated to the chiral limit are slightly different from those 
of ref. H. There are several reasons for the differences: i) in order to fix the scale and the 



strange quark mass we have used the lattice-plane method of ref. |25]; ii) in the present 
analysis, we use the "lattice dispersion relation" of eq. (|l8|), instead than the continuum one 
E'^ = m? ± X]j=i 3P? which was adopted in our previous study 0; iii) in order to reduce the 
systematic effects due to higher order terms in the chiral expansion, i.e. to higher powers of 
p-q, we have not used the results corresponding to p = 2tt/L{1, 0, 0) and q = 2it/L{—1, 0, 0). 
This choice stabilizes the results for Bi between /? = 6.0 and /3 = 6.2 whilst the results for 
the other ^-parameters remain essentially unchanged. 

We now describe the criteria followed in order to obtain our best estimates of the B- 
parameters. Although we have data at two different values of the lattice spacing, the statis- 
tical errors, and the uncertainties in the extraction of the matrix elements, are too large to 
enable any extrapolation to the continuum limit a — > : within the precision of our results 
we cannot detect the dependence of 5-parameters on a. For this reason, we estimate the 
central values by averaging the i?-parameters obtained with the physical mass at the 

two values of /3. Since the results at /? = 6.0 have smaller statistical errors but suffer from 



^ Although the continuum and lattice dispersion relations are equivalent at the order in a at which we are 
working, the latter gives a better description of the data for large momenta, see for example ref. [E6|. 
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Bi 


B2 


Bz 


Ba 


^5 


0.31 


0.67 ±0.21 


0.72 ±0.10 


0.90 ±0.40 


0.99 ±0.12 


0.21 ±0.24 


0.62 


0.68 lb 0.21 


0.63 ±0.06 


0.94 ±0.16 


0.98 ±0.08 


0.46 ±0.13 


0.96 


0.68 ± 0.21 


0.60 ± 0.04 


0.91 ±0.11 


0.95 ± 0.07 


0.54 ± 0.10 


1.27 


0.68 ±0.21 


0.58 ± 0.04 


0.89 ±0.08 


0.93 ± 0.07 


0.56 ±0.09 


1.38 


0.68 ±0.21 


0.58 ± 0.04 


0.88 ±0.08 


0.92 ±0.07 


0.55 ±0.09 


1.85 


0.67 ±0.21 


0.57 ±0.03 


0.84 ±0.07 


0.91 ±0.06 


0.56 ±0.08 


2.46 


0.68 ± 0.21 


0.56 ± 0.03 


0.82 ± 0.06 


0.91 ± 0.06 


0.58 ± 0.07 


4.00 


0.69 ± 0.21 


0.54 ± 0.03 


0.78 ± 0.05 


0.91 ± 0.06 


0.60 ± 0.07 



Table 3: Values of the B-parameters in the chiral limit, at different scales ji'^a^ for (3 = 6.2. 





Bi 


B2 


Bs 


^4 


B5 


0.31 


0.72 ± 0.14 


0.74 ± 0.06 


1.57 ±0.22 


1.19 ±0.06 


1.70 ±0.19 


0.62 


0.70 ±0.15 


0.69 ± 0.03 


1.22 ±0.10 


1.08 ±0.04 


0.92 ±0.09 


0.96 


0.70 ±0.15 


0.66 ± 0.03 


1.12 ±0.07 


1.05 ±0.03 


0.79 ± 0.06 


1.27 


0.70 ±0.15 


0.65 ±0.02 


1.06 ±0.06 


1.03 ±0.03 


0.79 ±0.05 


1.38 


0.70 ± 0.15 


0.64 ± 0.02 


1.06 ± 0.05 


1.02 ± 0.03 


0.79 ± 0.05 


1.85 


0.71 ±0.15 


0.63 ± 0.02 


1.02 ±0.05 


1.01 ±0.03 


0.77 ± 0.04 


2.46 


0.72 ±0.15 


0.61 ±0.02 


0.99 ± 0.04 


1.02 ±0.03 


0.77 ±0.04 


4.00 


0.74 ±0.15 


0.59 ±0.02 


0.90 ±0.03 


1.01 ±0.03 


0.81 ±0.03 



Table 4: Values of the B-parameters at the physical kaon mass, at different scales y^a^ for {3 — 6.0. 





Bi 


B2 


S3 


S4 


Bf> 


0.31 


0.67 ±0.21 


0.75 ± 0.07 


1.05 ±0.29 


1.05 ± 0.09 


0.57 ±0.18 


0.62 


0.68 ± 0.21 


0.66 ± 0.04 


0.98 ±0.12 


1.01 ± 0.06 


0.67 ±0.10 


0.96 


0.68 ± 0.21 


0.63 ± 0.03 


0.95 ± 0.08 


0.99 ± 0.06 


0.70 ± 0.08 


1.27 


0.68 ±0.21 


0.61 ±0.03 


0.92 ± 0.07 


0.97 ±0.05 


0.71 ±0.07 


1.38 


0.68 ±0.21 


0.61 ± 0.03 


0.91 ±0.06 


0.97 ±0.05 


0.70 ± 0.07 


1.85 


0.67 ±0.21 


0.60 ± 0.03 


0.88 ± 0.06 


0.97 ± 0.05 


0.70 ± 0.06 


2.46 


0.68 ± 0.21 


0.59 ± 0.03 


0.86 ± 0.05 


0.97 ± 0.05 


0.71 ± 0.06 


4.00 


0.69 ± 0.21 


0.57 ± 0.02 


0.82 ± 0.04 


0.97 ± 0.05 


0.73 ± 0.06 



Table 5: Values of the B-parameters at the physical kaon mass, at different scales if^a"^ for (3 — 6.2. 
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larger discretization effects, we do not weight the averages with the quoted statistical errors 
but take simply the sum of the two values divided by two. As far as the errors are concerned 
we take the largest of the two statistical errors. This is a rather conservative way of esti- 
mating the errors. In order to compare the results of Run A and Run B, we must choose 
the same physical renormalization scale fi. Using the estimates of the lattice spacing given 
in table |l], we have taken /i^a^ = 0.96 and fj.'^a'^ = 0.62, corresponding to /i = 2.08 GeV and 
fi = 2.12 GeV, at /? = 6.0 and 6.2 respectively. We quote the results as obtained at /U = 2 
GeV, since the running of the matrix elements between ^ 2.1 and 2.0 is totally negligible 
in comparison with the final errors. In table we summarize the values which have been 
used to give the final estimates. The columns denoted by have been used to get the 

final results in eq. (|3|). 



-B(/i ~ 2 GeV) 


rriK = 


rriK = 


exn 

rriK = mx 


niK = 


exp 

rriK = 




13 = 6.0 


P = 6.0 


(3 = 6.0 


/3 = 6.2 


(3 = 6.2 




ref.i 


this work 


this work 


this work 


this work 


Bi 


0.66(11) 


0.70(15) 


0.70(15) 


0.68(21) 


0.68(21) 


B2 




0.61(3) 


0.66(3) 


0.63(6) 


0.66(4) 


Bs 




1.10(8) 


1.12(7) 


0.94(16) 


0.98(12) 


Bi 


1.03(3) 


1.04(4) 


1.05(3) 


0.98 (8) 


1.01(6) 


B5 


0.72(5) 


0.68(7) 


0.79(6) 


0.46(13) 


0.67(10) 



Table 6: B-parameters at the renormalization scale = a ""^ ^ 2 GeV, corresponding to /i^a^ 
0.96 and fi'^a'^ = 0.62 at (3 = 6.0 and 6.2 respectively. All results are in the RI (MOM) scheme. 



In ref. B2 and i^a have been obtained at (3 = 6.0 with the Wilson action and the 
operators renormalized perturbatively in the MS scheme; the result is 



B2 = 0.59 ±0.01 
B3 = 0.79 ±0.01 



(22) 



Although a direct comparison is not possible (our results are in the RI scheme), to the extent 
that the matching coefficients between the two schemes are a small effect [^, comparison of 
eqs.(^ and (22) suggests that perturbative renormalization gives significantly different results 
in some cases. This confirms the need for non perturbative renormalization. 



5 Conclusions 

In this paper we have presented a lattice calculation of the matrix elements of the most general 
set of AS* = 2 dimension-six four-fermion operators, renormalized non-perturbatively in the 
RI (MOM) scheme @]-0- The calculations have been performed at two different values of 
the lattice spacing a. Although our precision is not sufficient to make an extrapolation to the 
continuum limit, the comparison between the results on two different lattices allows a better 
estimate of the final errors. The main results for the five S-parameters are summarized in 
tab. 1^. From this table, we have extracted our best estimates which are given in eq. (^). 
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We observe that the lattice values of 4 are close to their VSA whereas this is not true 
for -Bi,2,5- Our results allow an improvement in the accuracy of phenomenological analyses 
intended to put bounds on basic parameters of theories beyond the Standard Model. 
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